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Abstract. 

The precise knowledge of the atmospheric neutrino fluxes is a key ingredient 
in the interpretation of the results from any atmospheric neutrino experiment. In 
the standard atmospheric neutrino data analysis, these fluxes are theoretical inputs 
obtained from sophisticated numerical calculations. In this contribution we present 
an alternative approach to the determination of the atmospheric neutrino fluxes based 
on the direct extraction from the experimental data on neutrino event rates. The 
extraction is achieved by means of a combination of artificial neural networks as 
interpolants and Monte Carlo methods. 
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1. Introduction 

One of the most important breakthroughs in particle physics, and the only solid evidence 
for physics beyond the Standard Model, is the discovery - following a variety of 
independent experiments pQ - that neutrinos are massive and consequently can oscillate 
among their different flavor eigenstates. The flavor oscillation hypothesis has been 
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supported by an impressive wealth of experimental data, one of the most important 
pieces of evidence coming from atmospheric neutrinos. 

Atmospheric neutrinos are originated by the collisions of cosmic rays with air nuclei 
in the Earth atmosphere. These neutrinos are observed in underground experiments 
using different techniques 0. In particular in the last ten years, high precision and 
large statistics data has been available from the SuperKamiokande experiment which 
has clearly established the existence of a deficit in the //-like atmospheric events with 
the expected distance and energy dependence from — > v T oscillations with oscillation 
parameters Am 2 ~ 2 x 10~ 3 eV 2 and tan 2 9 = 1. This evidence has also been confirmed 
by other experiments, both athmosferic such as MACRO and Soudan 2 and accelerator- 
based like MINOS and K2K. 

The expected number of atmospheric neutrino events depends on a variety of 
components: the atmospheric neutrino fluxes, the neutrino oscillation parameters and 
the neutrino-nucleus interaction cross section. Since the main focus of atmospheric 
neutrino data interpretation has been the determination of neutrino oscillation 
parameters, in the standard analysis the remaining components of the event rate 
computation are inputs taken from other sources. In particular, the fluxes of atmospheric 
neutrinos are taken from the results of numerical calculations, like those of Refs. 011], 
which make use of a mixture of primary cosmic ray spectrum measurements, models for 
the hadronic interactions and simulation of particle propagation 0. 

The attainable accuracy in the independent determination of the relevant neutrino 
oscillation parameters from non-atmospheric neutrino experiments makes it possible to 
attempt an inversion of the strategy in the atmospheric neutrino analysis: to use the 
oscillation parameters as inputs in the data analysis in order to extract from data the 
atmospheric neutrino fluxes. 

There are several motivations for such direct determination of the atmospheric 
neutrino flux from experimental data. First of all it would provide a cross-check of 
the standard flux calculations as well as of the size of the associated uncertainties 
(which being mostly theoretical are difficult to quantify). Second, a precise knowledge 
of atmospheric neutrino flux is of paramount importance for high energy neutrino 
telescopes , both because they are the main background and they are used for detector 
calibration. Finally, such program will quantitatively expand the physics potential 
of future atmospheric neutrino experiments. Technically, however, this program is 
challenged by the absence of a generic parametrization of the energy and angular 
functional dependence of the fluxes which is valid in all the range of energies where 
there is available neutrino data. 

In this contribution we present the first results on this alternative approach to the 
determination of the atmospheric neutrino fluxes: we will determine these fluxes from 
experimental data on atmospheric neutrino event rates, using all available information 
on neutrino oscillation parameters and cross-sections. The problem of the unknown 
functional form for the neutrino flux is bypassed by the use of neural networks as 
interpolants, since they allow us to parametrize the atmospheric neutrino flux without 
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having to assume any functional behavior. Additional details on our approach can be 
found in jjj. 

Indeed the problem of the deconvolution of the atmospheric flux from experimental 
data on event rates is rather close in spirit to the determination of parton distribution 
functions in deep-inelastic scattering from experimentally measured structure functions 
jH]. For this reason, in this work we will apply for the determination of the atmospheric 
neutrino fluxes a general strategy originally designed to extract parton distributions in 
an unbiased way with faithful estimation of the uncertainties^ [HI H21 HHj • 

2. General strategy 

The general strategy that will be used to determine the atmospheric neutrino fluxes was 
first presented in Ref. [TT] (see also [Ej). It involves two distinct stages in order to go 
from the data to the flux parametrization. In the first stage, a Monte Carlo sample of 
replicas of the experimental data on neutrino event rates ( "artificial data" ) is generated. 
These can be viewed as a sampling of the probability measure on the space of physical 
observables at the discrete points where data exist. In the second stage one uses neural 
networks to interpolate between points where data exist. In the present case, this second 
stage in turn consists of two sub-steps: the determination of the atmospheric event rates 
from the atmospheric flux in a fast and efficient way, and the comparison of the event 
rates thus computed to the data in order to tune the best-fit form of input neural flux 
distribution ( "training of the neural network" ) . 

In the present analysis we use the latest data on atmospheric neutrino event rates 
from the Super Kamiokande Collaboration [2]. Higher energy data from neutrino 
telescopes like Amanda pH] is not publicly available in a format which allows for its 
inclusion in the present analysis and its treatment is left for future work. 

The latest Super Kamiokande atmospheric neutrino data sample is divided in 9 
different types of events: contained events in three energy ranges, Sub-GeV, Mid- 
GeV and Multi-GeV electron- and muon-like, partially contained muon-like events and 
upgoing stopping and througoing muon events. Each of the above types of events is 
divided in 10 bins in the final state lepton zenith angle fa. Therefore we have a total of 
-^dat = 90 experimental data points, which we label as 

Rt P \ z = l,...,iV dat . (1) 

Note that each type of atmospheric neutrino event rate is sensitive to a different region 
of the neutrino energy spectrum. In particular we note that for energies larger than 
E v ~ few TeV and smaller than E v ~ 0.1 GeV there is essentially no information on 
the atmospheric flux from the available event rates. 

The purpose of the artificial data generation is to produce a Monte Carlo set of 
'pseudo-data', i.e. iV re p replicas of the original set of iV dat data points, Rr 1 *'^ such 

| This strategy has also been successfully applied with different motivations in other contexts like tau 
lepton decays jH] and B meson physics ^0] 
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that the N rep sets of N dat points are distributed according to an Ad^-dimensional 
multi-gaussian distribution around the original points, with expectation values equal 
to the central experimental values, and error and covariance equal to the corresponding 
experimental quantities. This is achieved by defining 



R 



(art) (A:) 



o(exp) + Jk) tot 
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where N Tcp is the number of generated replicas of the experimental data, and where 
r\ are univariate gaussian random numbers with the same correlation matrix as 
experimental data. 

In our case each neural network parametrizes a neutrino flux. Due to available 
experimental precision, in this work we will assume the zenith and type dependence of 
the flux to be known with some precision and extract from the data only its energy 
dependence. That is, if NN (E„) is the neural network output when the input is the 
neutrino energy E u , then the neural flux parametrization will be 

^ (E„, c„, t) = NN (E v ) $( rcf ) {E v , c, t) . (3) 

is a reference differential flux, which we take to be the most recent computations 
of either the Honda [Sj or the Bartol jlj collaborations, which have been extended to 
cover also the high-energy region by consistent matching with the Volkova fluxes. 

Now let us describe a fast and efficient technique to evaluate the atmospheric 
neutrino event rates i?| net ^ for an arbitrary input neural network atmospheric flux 
,t). As shown in [2j it is possible to write the theoretical predictions as 



a sum of the elements of a bin-integrated flux table \l/ 
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where the coefficients C\ ztl which are the most time-consuming ingredient, need 
only to be precomputed once before the training, since they do not depend on the 
parametrization of the atmospheric neutrino flux. 

The determination of the parameters that define the neural network, its weights, 
is performed by maximum likelihood. This procedure, the so-called neural network 
training, proceeds by minimizing an error function, which coincides with the x 2 °f the 
experimental points when compared to their theoretical determination obtained using 
the given set of fluxes: 
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The x 2 has to be minimized with respect to {uji}, the parameters of the neural network. 
The minimization of Eq. (0) is performed with the use of genetic algorithms ^3] . 

Thus at the end of the procedure, we end up with A^p fluxes, with each flux 

$(net)( fc ) 

given by a neural net. The set of A^p fluxes provide our best representation of 
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Figure 1. Comparison of the 1 — a uncertainty band of our parametrization of the 
atmospheric neutrino flux with different atmosferic neutrino flux computations and 
with the latest AMANDA data. 



the corresponding probability density in the space of atmospheric neutrino fluxes: for 
example, the mean value of the flux at a given value of E v is found by averaging over 
the replicas, and the uncertainty on this value is the variance of the values given by the 
replicas. 

3. Results and discussion 

Now we discuss the results of our determination of the atmosferic neutrino fluxes. Note 
that any functional of the neutrino flux (like for example event rates and its associated 
uncertainty) can be computed taking the appropriate moments over the constructed 
probability measure 

iVrcp 

(^($(^))) rep = T^E^( <|,(nCt)(fc) (^)) • ( 6 ) 

re p k=l 

In Fig. ^ we show our results for the atmosferic neutrino flux as compared with the 
Honda and Bartol computations as well as with some direct measurements from Amanda 
[To] . Note that above 1 TeV we are in the extrapolation region. We show in Fig. |2]how 
our neural network parametrization successfully reproduces the features of experimental 
data. A more detailed discussion of the results can be found in [7] 

The work presented in this contribution can be extended in several directions. 
First of all one could reduce the uncertainty in the large E v region by incorporating 
in the fit atmosferic neutrino measurements from neutrino detectors like Amanda 
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Figure 2. Comparison of the neural network results with experimental data for the 
Multi-GeV data sample. 



[IBJ HH] . Second, we could use at even larger energies simple parametrizations of the 
flux to obtain a better estimation of backgrounds at neutrino telescopes. Finally, we 
could estimate the required statistics that would be required in forthcoming atmosferic 
neutrino experiments in order to determine from data also the zenith angle and flavor 
dependence of the neutrino flux. 
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